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Abstract 

We study thermal fluctuations and capacitive effects on small Joseph- 
son Junction Rings (JJR) that mimics the rectification phenomena recently 
observed in triangle shaped mesoscopic superconductors, due to the superpo- 
sition of the field induced persistent current with the bias current. At finite 
temperature we predicted that the amplitude of the rectified signal depends 
strongly on the current contacts configuration on the JJR, in coincidence 
with experiments. In addition we analize the range of parameters where a 
closed loop of capacitive junctions is an appropriate model to explain the 
experimental observations. We conclude that the closed loop of weak links, 
a JJR, is a simple, robust and good enough model to explain the observed 
voltage rectification effects on mesoscopic superconducting samples for a wide 
range of temperature. 
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1. INTRODUCTION 



Technological progress allowing the manipulation of ratchet effects in 
complex systems such as cold atoms [JJ, colloids [2J, cells [3J, fluids [I], elec- 
trons in semiconductors [5], extended magnetic walls [6j, vortices in Joseph- 
son systems [3, E], motivate this work. In all the mentioned examples the 
flux of driven particles interact with an asymmetric potential. The key ques- 
tion on this field is that it is not simple in each case to know feasibly which 
are the mechanisms underlying these many bodies rectification effects. It 
is well known that the velocity of a single overdamped particle under an 
applied ac force can be rectified by introducing a periodic potential with 
broken reflection symmetry [9j. We can thus ask whether this effect, known 
as rocking ratchet (because the system is rocked with an external ac force), 
can be still observed as a rectified dc voltage response in transport mea- 
surements in singly connected mesoscopic superconductors under external 
ac bias currents, thermal fluctuations and external homogenous magnetic 
field applied. Indeed finite dc voltage was measured under zero averaged 
ac current applied and it was extensively studied in superconducting (SC) 
nanostructures [TQl [TT], [12] . But this system differs from examples mentioned 
previously, because the essential broken symmetry required to obtain a recti- 
fier is controlled differently in these samples, it is controlled by its geometry 
(sample shape) and the configuration of contacts leads (position of the cur- 
rent/voltage probes). In particular, motivated by these recent experiments 
reporting original voltage rectification effects in mesoscopic superconducting 
triangles [TT] we simulated a more realistic model for such system in order 
to understand and corroborated the physics behind the observed ratchet ef- 
fects, with special emphasis in the finite temperature effects not taking into 
account in previous and simpler models. 

In a pionner analytical work [13] different asymmetric type-II supercon- 
ducting structures (rings and strips) were already proposed and analyzed 
under current and magnetic field applied. They show how magnetic flux 
penetration depending on the geometric parameters is fieldlike or current- 
like. Later on, measurements were realized in asymmetric SC rings [TJJ IT5] . 
In addition similar rectification effects were reported in a singly connected 
structured if the current injection is off-center [10] . From mentioned works it 
is widely accepted that the effects observed are due to the asymmetry caus- 
ing a difference in critical current for a positive or negative applied external 
current which is compensated or reinforced by the field induced persistent 
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current. Further an important parameter, temperature, should be taking into 
account to model in detail the phenomena observed. First because all results 
reported in single triangle samples are done close to T c [101 [UJ [12] , it means 
where thermal fluctuations are important and far away from T = where 
previous theoretical analysis was done [TT] . In addition it is well known that 
geometry in microsized or in nanosized samples could play essential roles in 
fluctuations phenomena and even more, very recently it was shown that the 
superconducting order parameter present larger thermal fluctuations in the 
corners of sharp shaped samples |16j . At present, when smaller dimensions 
are manipulated it is motivating to study effects of thermal fluctuations even 
at temperature much lower than T c where dimensionality start to compete 
and the limits between thermal fluctuations and quantum fluctuations are 
not clear. In short, to understand how thermal fluctuations influence the 
rectification phenomena at all finite temperature is important. 

In this paper we model voltage rectification effects observed in thin meso- 
scopic superconducting triangles [TT] using two kind of small Josephson Junc- 
tion Rings: a resistive one, using the RSJ (resistively shunted junction) model 
and a capacitive model, using the RCSJ (resistevely and capacitive shunted 
junction) model. We introduce finite temperature effects adding white noise 
to the resulting dynamics equations within the RSJ model. We find good 
qualitative agreement with experiments. 

2. SIMULATED MODEL 

We study at finite temperature the dynamics of small Josephson Junction 
Rings (closed loops) with a finite number of identical Josephson junctions 
(JJ): three J J N = 3, four JJ N = 4 or five J J, N = 5 as sketched in 
Figjl] for different positions of the current injection. An external sinusoidal 
ac current is applied in different loop points as shown in Figjl] in order to 
mimics the different contact configurations, samples A, B and C used in 
experiments [TT]. 
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Sample A Sample B Sample C 




Figure 1: (Color online) Schematic small Josephson junctions rings modeling 
superconducting triangles, using either N = 3 JJ, N = 5 JJ or iV = 4 JJ 
(upper panels), according to the experimental samples A, B or C with differ- 
ent contact positions (lower panels). Note the different asymmetric current 
injection for N = 3, A, below the geometrical center and N = 5, B, above the 
geometrical center in opposite to the symmetric current injection for N = 4, 
C. Lower panels (generosity of Nele Schildermans and A.V. Silhanek) are 
scannig electron microscopy images of the superconducting Al triangles used 
in experiments. 

In Fig|2]we show a schema of the shape in which superconductivity nu- 
cleation occurs in triangle shaped mesoscopic samples. It is well known that 
nucleation in mesoscopic samples is dominated by the surface superconduc- 
tivity and its shape and boundaries, i. e by its geometry. In fact, it was clearly 
shown that surface nucleation is enhanced in wedge shaped superconductors 
[TTI H~8] . then the idea to model the triangle as a loop of strong supercon- 
ducting electrodes (SC islands) at the triangle corners connected by weak 
links is reasonable. Those weak links, could be an insulator as Josephson 
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proposed originally (SIS J J) or a normal metal (SNS JJ), or simpler a short 
narrow constriction of the same superconductor (ScS JJ) [20]. Our main 
assumption for the calculations is that the superconducting order parameter 
p(r) = \p(r)\e ( T > is such that |p(r)| = p with p the same constant on all 
SC islands and 8(r) is spatially constant in each island. Then our variables 
will be the superconducting phases in each SC island. 



K A 




T > T c T ~T C 



Figure 2: (Color online) At large tempererature the Al triangles are normal 
samples but at T ~ T c nucleation of the superconducting condensate starts 
from the sharp sample corners. The order parameter \l/ is maximum at the 
triangle vertices and minimum at the middle of the sides, then a ring of weak 
links (JJ) is used to model this system. 

We started solving the simpler case, the ring of weak links of SNS junc- 
tions, using the resistively shunted junction (RSJ) model plus thermal fluc- 
tuations and solving numerically the superconducting phase dynamics using 
Langevin dynamics technique [19] . We considere a total magnetic field B 
applied perpendicular to the ring, spatially and temporally constant. The 
Hamiltonian of a closed loop of N SNS-j unctions is the following: 

N-l 

H = -Ej ^ cos (0n - O n ) (1) 
n=0 

where Ej is the Josephson energy, n = 6*(r n )— 6*(r n _i) is the superconducting 
phase difference at the junction n = 0, ...,N, and 8{r n ) is the phase of the 
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superconducting island centered at 



rn = R (- cos ( x + sin ( ^ j y) (2) 

being R the ring radius. Note here that our loops imitating the triangle sam- 
ples, are ideal junction circles where the needed triangle asymmetry can be 
introduce controlling the number of junctions and the points where external 
currents are injected and extracted (see Fig. Q. The total magnetic field is 
B = V x A. The magnetic field contribution to the phase difference a n is 
the line integral of the vector potential between sites n and n — 1, 
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A(l).dl, (3) 



being $ = h/2e the single superconducting quantum flux. After taking 
the magnetic field applied perpendicular to the sample plane, B = Bz, the 
invariant gauge A = Bxy for the vector potential A and the flux quantum 
number through the ring $/$o — TtR 2 B/§Q we obtain after solving the 
corresponding previous integral in polar coordinates: 



2$ r 
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(4) 



The current flowing in the JJ between two superconducting islands is 
modeled by a typical electric circuit build up of an ideal JJ in parallel with 
a normal resistance, i.e by the sum of the Josephson supercurrent and the 
normal electron current. Conservation of current is considered with Kirchhoff 
laws fulfilled in each node. We inject a current I between junctions N — 1 
and 0, and extract it S junctions away, between junctions 5 — 1 and 5. The 
resulting set of dimensionless equations for the currents flowing in the ring 
is the following: 

<fin = hp - sin(0„ - a n ) + r(n, t), < n < 5 - 1 (5) 
4> n = hp - I - sin(0 n - On) +T(n,t), 5 < n < N - 1 (6) 

1 N ^ 

I up = (1 - S/N)I + sin ^« _ a «) + r ( n ' *)■ ( ? ) 

n=0 
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which are N first order differential equations for the time evolution of the 
N phase variables {4> n }n=o- ^"p * s the current flowing in the upper branch 
of the ring, up with respect to the injection and extraction current points, 
and I — I up = I down- Let us note that each junction interacts with all the 
others through I U p({4>n}n=o) > the total current in the upper branch of the 
circuit, which represents a kind of mean-field interaction plus a drive. In 
addition finite temperature effects are taking into account through an added 
Langevin noise term T, introduced as current fluctuations, which models the 
contact with a thermal bath at temperature T and satisfies the condition for 
non-correlated white noise: 

Ik T 

(T(n, t)T(n', t')) = -^5 n , n ,5(t - t'). (8) 

where is the normal state resistance. In short the model control 
parameters are the external current J, magnetic field applied or magnetic 
flux $/5>o and temperature T. In addition one can play with the source- 
drain distance 5 and the number of junctions N. 

Langevin dynamical equations Eq. |5|6] using conditions Q and (IT]) are 
solved numerically using a second order Runge-Kutta-Helfand-Greenside al- 
gorithm under an external sinusoidal current applied I = I ac sin(27iut). We 
calculate the mean voltage response computing the instantaneous dc voltage 
drop v between source and drain: 

5-1 5-1 

v = ^2<p n = SI up - ^ sin (0n - a n ) + r(n, t). (9) 

n=0 rc=0 

In fact there is a second equivalent form to compute the voltage drop, if 
upper and lower paths contain a different number of junctions: 

AT-l N-1 

v = -J2<l>n = {N-d)(I- I up ) + sin(0„ - a n ) + F(n,t). (10) 

n=5 n=5 

We normalize currents by the single junction critical current Jo, voltages 
by -Rat Jo > temperature by Ej/ks and the unit of time is Tj=2iicRiyI /^ . 

Possibly weak links sketched in Fig. [2] could be more complex than the 
simpler SNS junctions, like tunnel superconducting-insulator-superconducting 
(SIS) junctions and a more complete description is required. Then we probe 



7 



if the essential features of previous experiments [HU ED] are as well described 
by a RCSJ (resistively and capacitively shunted junction, [20]) ring. In this 
model the weak link ( JJ) is modeled by an ideal JJ shunted in parallel by a re- 
sistance R c and a capacitance C. The addition of capacitive effects accounts 
for the geometric shunting capacitance between the two superconducting is- 
lands. Now the following set of second order differential equations should be 
solved: 

4> n + v4>n + sin((f) n - a n ) = I up , < n < 5 - 1 (11) 

4>n + V<fin + sin(0 n - a n ) = I up - I, 5 < n < N - 1 (12) 

being rj the main parameter now, a measurement of the system dissipa- 
tion, r) — 1/ a//3 and /3 is known as the McCumber parameter [21], or damping 
parameter equal to (u p R c C) 2 with u p the plasma frequency of the junction 

= y/2eI /ftC. 

3. RESULTS 

In order to study the influence of contacts on the rectification effects 
observed in superconducting triangles [TT] the model for RSJ rings described 
in previous section was first solved numerically. We calculated the rectified 
mean dc voltage Vd c = (v) as a function of the ac current amplitude, I ac , in the 
low frequency limit at finite temperatures and for different magnetic fields 
applied perpendicularly to the samples. The results for closed loops with 
N = 3 junctions (which mimics the experimental contacts configuration of 
sample A) and N = 5 (sample B), both with the same source-drain distance 
5 = 2, are shown in Fig. |3j 

We can clearly see that at finite temperature both, the N = 3 and N = 5 
devices display finite rectified voltage, i.e. \Vd c \ > above a critical cur- 
rent. Examples for a fixed magnetic field corresponding to n = —1, i.e. for 
^/c&o — —1/4 and two temperature values are shown. The junction loop with 
N = 5 presents a clear difference comparing with the case N = 3, the rectified 
voltage is inverted, i.e. Vd c < 0, resulting both junction loops in good choices 
to reproduce the experimental results obtained in Ref. [11] for both type of 
asymmetric contact configurations, sample A and sample B, for which in- 
verted rectified voltage responses were reported. In short, this simple model 
including thermal fluctuations mimics correctly recent experimental results 
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Figure 3: (Color online) Rectified voltage vs ac-current amplitude applied for 
N = 3 (sample A) and T = 5 x 10~ 7 (black ■), iV = 3 and T = 6 x l(T 6 (red 
•),N = 5 (sample B) and T = 5 x 10~ 7 (green ▲), N = 5 and T = 6 x 1(T 6 
(blue -jc). For all examples shown: $/$o = —1/4. The symmetric sample C, 
N — 4, it is not shown for clarity because Vd c = for all I ac and all magnetic 
field applied. 

[TT] . Other details could be extracted from Fig.|3j all examples shown present 
a critical current amplitude threshold that decreases with temperature and 
it is larger for the device with the larger number of junctions, iV = 5. We 
can also conclude that the maximum of | Vd c \ is almost the same in both cases 
when the temperature is closer to zero, but when temperature is increased 
these maximum decrease, as it is expected when thermal fluctuations are 
present, but the rates are different depending on the sample type, the max- 
imum rectified voltage |Vd c | maa; decreases quicker in the larger loop, N = 5. 
At T = 5 x 1CT 7 both cases present a maximum at |Vrf c | maa; ~ 0.21 but when 
temperature is decreased more than an order of magnitude, T = 6 x 10 -6 , 
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the maximum is |V^ c | max ~ 0.15 for N = 3 and |Vrf c | maa; ~ 0.1 for N = 5. 
In addition \Vd c \ decays to zero with I ac , as expected, but this behavior is 
quicker for N = 5. 
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Figure 4: (Color online) Oscillations of the rectified dc-voltage V^ c vs mag- 
netic flux $/$o, comparing both samples and two finite temperatures: N = 3 
and T = 5 x 10~ 7 (black ■), N = 3 and T = 1 x 10~ 5 (red •), JV = 5 and 
T = 5 x lO" 7 (blue *), N = 5 and T = 1 x 10~ 5 (green a). 

Another interesting feature to analize is the voltage response related to 
the total magnetic field applied, B. An important behavior is obtained and 
analyzed in Figj4j, where qualitatively the same oscillations for V^ c as a func- 
tion of vorticity (magnetic flux $/$o through the ring) are also obtained in 
the experiments [11]. Three characteristics of the oscillations dependency on 
magnetic field are found to be independent of temperature: (a) the recti- 
fication effects disappear if $/$o = ra /2, being n an integer, (b) the max- 
imum absolute value of the rectified voltage appears always the condition 
$ = (2ra + l)$ /4 is fulfilled and (c) \V dc \ ^ if $/$ ^ n/2. The observed 
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Figure 5: (Color online) Rectified dc voltage V^ c vs temperature T in log- 
normal scale for I ac = 2 in all examples and N = 3 and $/$o = —1/4 (red 
•), iV = 3 and $/$ = -1/8 (black ■), N = 3 and $/$ = -1/2 (blue *), 
iV = 5 and $/$ = -1/4 (green A). 

dependences of rectified voltage on magnetic field are in accordance with ex- 
pected Little-Parks oscillations, which occur every time the vorticity of the 
system increase due to the entrance of a new vortex: L — > L + 1 [22J. In ad- 
dition a noticeable feature appears, the sign of the rectified voltage changes, 
from negative if n < $/$o < n + 1/2 to positive if n+ 1/2 < $/$o < n + 1 
(being n an integer) for N = 3 and signs are inverted for N = 5. This oppo- 
site behavior between both devices is related with its asymmetry and it is the 
best evidence that indeed the observed rectification effects are due to the su- 
perposition of the external applied sinusoidal current (bias current) with the 
field induced persistent current circulating around the JJRs. The external 
current injection is well below the geometrical center for iV = 3 device on the 
contrary to iV = 5 device, then an opposite current compensation behavior 
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was obtained. It is worth to mention that all simulations were performed as 
well for the symmetric current injection case, sample C, and no voltage re- 
sponse was observed, corroborating the essential need of an off-center current 
injection (system asymmetry induced externally for the contacts configura- 
tion). Note that \Va c \ decreases to zero while temperature is increased, in 
both cases. In order to analyze this last feature in detail comparing both 
contact configurations, the rectified dc voltage as a function of temperature 
is shown in Fig. [5] for I ac = 2 and some typical values of the magnetic field 
applied, corresponding to $/$o — (2n + l)/4 = —1/4 (maximum voltage 
response), $/$ = (4ra + l)/8 = —1/8 (intermediate voltage response) and 
$/$o — — 1/2 (minimum voltage response). For fields with finite response, 
the absolute value of the rectified voltage decrease smoothly to zero with 
similar temperature functionality while temperature is incremented slowly in 
several orders of magnitude. It is worth to mention that it remains finite 
for a wide range of temperature. Analyzing the three curves presented for 
iV = 3, a rectified voltage tendency to zero as the magnetic flux approach to 
n/2 is observed as well, in coincidence with previous Fig. |4| These results 
show clearly the robustness of the rectification effects under a wide range of 
thermal fluctuations for all magnetic field applied such that $/$o 7^ n/2. 

In short we found that simple loop of weak links modeled with a resis- 
tive shunted model (RSJ) plus thermal fluctuations is in good agreement 
with experimental results. The following question naturally arises: is a ring 
of capacitive Josephson junctions, i.e. a capacitive and resistive shunted 
model (RCSJ) a suitable model to explain the observed phenomena?. In 
order to get an insight into this point we include the capacitive effects in 
our simulations. As an illustration in Figj6] a comparison of the mean dc 
voltage response V& c vs I ac among the three type of contact configurations 
(N = 3,4, 5.) is shown. The results presented were obtained at zero temper- 
ature for a magnetic flux $/$o = (2n + l)/4. A capacitive model without 
thermal fluctuations is shown to be a good aproximation to the experimental 
observations if /3 < 1. In Fig. [6] one /3-value is shown for clarity, (3 = 0.01, 
but simulations were performed incrementing j3 till unity, and similar results 
were obtained regarding our previous main conclusions in comparison with 
experimental results observed in Ref. jllj . 
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Figure 6: Rectified voltage response obtained using RCSJ model, as a func- 
tion of ac current amplitude for (5 = 0.01: N = 3 (red ■); N = 4 (green •); 
N = 5 (blue *). For all examples $/$ = -1/4. 

4. CONCLUSIONS 

In summary, at finite temperature we predicted within a model system 
that the amplitude of the rectified signal depends strongly on the current 
contacts configuration in the JJR, in good agreement with recent diode ef- 
fects measured in Al mesoscopic triangles [TQl EH] and in coincidence as well 
with a more sophisticated model including capacitive effects. Our results 
corroborate that the clue behind the observed rectification effects in triangle 
shaped SC is indeed the superposition of the field induced persistent current 
with the bias current, being these phenomena robust and strong enough in 
front of thermal fluctuations and for more complicated modeling of the weak 
links. The key ingredients to observe the dc voltage rectification described 
here under thermal fluctuations are indeed field induced persistent currents 
and asymmetry induced by an off-center injection of external currents. This 
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recipe suggests that similar rectification effects, strongly robust under finite 
temperature, should be also present in every system with a persistent current 
and an asymmetric current path, being useful to design and control newer 
and smaller diode devices. 
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